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ABSTRACT 


The  coulomb  repulsion  between  the  electrons  Is 
represented  as  the  sum  of  a  long-range  unperturbed 
part  and  a  short-range  perturbation.   Wang's  form 
of  the  Heltler-London  approximation  is  used  as  the 
approximate  eigenfunction  for  the  unperturbed  hamll- 
tonian.   The  correction  for  the  short-range  per- 
turbation Is  made  by  a  method  based  on  the  Impulse 
approximation . 


Present  address,  Stevens  Institute  of  Technology, 
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INTRODUCTION 

At  points  where  two  or  more  electrons  coincide,  the 
analytic  behavior  of  trial  functions  commonly  used  In 
Raylelgh-Rltz  calculations  differ  from  that  of  solutions 
of  Schrodlnger ' s  equation.   Consequently,  these  trial  func- 
tions are  poor  approximations  to  the  correct  elgenfunctlons 
In  domains  where  the  distance  between  two  electrons  Is 
sufficiently  small.   We  attempt  to  study  the  resulting  error 
In  the  Raylelgh-Ritz  energy  estimate  by  Improving  the  trial 
function  In  these  domains.   We  do  not  expect  to  achieve  this 
by  superposing  functions  with  the  incorrect  analytic  behavior, 
for  this  leads  to  slow  convergence.    '  Instead,  we  try  to 

satisfy  the  singularity  condition,  at  least  formally.   We 
separate  the  coulomb  repulsion  Into  two  parts,  a  long-range 
unperturbed  part  and  a  short-range  perturbation.   This  sepa- 
ration Is  arbitrary  and  contains  parameters  which  are  adjusted 
at  a  later  stage  of  the  calculation.   We  assume  the  ground 
state  elgenfunctlon  of  the  unperturbed  hamlltonlan  Is  approx- 
imated well  by  a  known  function,  such  as  the  Heltler-London, 
molecular  orbital,  or  small  superpositions  of  these.   We 
account  of  the  perturbation  by  an  adaptation  of  the  Impulse 
approximation,  ^  ~    which  applies  to  short-range  forces 
which  need  not  be  weak.   The  parameters  used  In  decomposing 
the  coulomb  repulsion  are  finally  adjusted  by  a  self-consistency 
condition.   The  formal  details  are  described  in  Section  1. 
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The  Integrals  needed  for  a  direct  variational  calculation 
are  difficult  to  evaluate.   To  circumvent  this  problem,  we  use 
a  simple  non-varlatlonal  condition  which  makes  the  ground 
state  energy  of  the  perturbed  and  unperturbed  hamlltonlans 
approximately  equal,  assiomlng  the  basic  model  Is  valid.   If 
this  condition  can  be  applied  to  adjusting  the  parameters  In 
the  unperturbed  hamlltonlan,  the  effect  of  the  perturbation 
would  be  small  and  the  difficult  Integrals  could  be  estimated 
crudely,  since  they  do  not  occur  in  the  zero-order  problem. 
The  condition  Involves  an  Integral  which  Is  much  simpler  than 
those  needed  for  the  variational  estimate,  but  the  evaluation 
still  requires  further  approximation.   Here  we  use  an  approx- 
imation which  simplifies  the  calculation  enormously,  and  also 
has  an  attractive  physical  interpretation.   However,  it  seems 
that  this  approximation  is  too  drastic.   Additional  study  of 
the  relevant  integral  seems  practical  as  well  as  useful. 
The  details  are  described  in  Section  11. 

To  have  a  definite  problem,  we  calculate  the  ground  state 
energy  of  the  hydrogen  molecule.   We  expect  the  unperturbed 
problem  to  be  relatively  simple,  while  the  effects  of  the  short- 
range  perturbation  probably  do  not  depend  on  the  details  of 
the  molecule,  but  only  on  the  general  features  of  the  un- 
perturbed eigenfunction.   Wang's  modification  of  the  Heitler- 

(7) 
London  trial  function    is  assumed  adequate  for  the  un- 
perturbed problem.   We  note  that  the  extension  to  problems 


Involving  more  than  two  electrons  is  straightforward,  similar 
to  the  generalization  of  the  Impulse  approximation.     This 
is  the  motivation  for  this  study  since  it  is  known  that 
electronic  energies  of  the  hydrogen  molecule  may  be  calculated 
with  high  accuracy  if  the  distance  between  electrons  is  used 
directly  in  the  trial  function.    ' 
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I.   GENERAL  THEORY 

If  the  protons  are  Infinitely  heavy,  the  hamlltonian  for 
the  hydrogen  molecule  is 


(1.1)        H  =  K  +  V^  +  i 


Here  K  is  the  kinetic  energy,   V„   the  attraction  between 
electrons  and  protons,   r   the  distance  between  electrons,  and 
atomic  units  are  used.   The  unperturbed  hamlltonian  is  obtained 
by  replacing   (l/r)   by  a  modified  repulsion  V   (r) 


(1.2)         H^  =  K  +  V^  +  V^^(r) 


(1.3)        V^^(r)  ^l-Z-f^+  ^   e-^^ 


The  parameters  A  and  [i     are  determined  later.   The  functional 
form  (1.3)  is  used  for  convenience,  and  eliminates  the  diver- 
gence in  the  energy  but  not  in  all  its  gradients.   This  could 

clearly  be  done  by  using  a  different  form  of  V   ,  namely  an 

2 
integral  function  of  r   rather  than  r.   Clearly  the  term  in 

(1.3)  proportional  to  [i      compensates  the  repulsion  removed  at 

small  r.   For  each  A  ,  we  expect  some  positive  m.  will  make 

the  ground  state  eigenvalues  of  H   and  H  equal.   We  now 

use  Wang's  form  for  the  spatial  part  of  the  unperturbed  trial 

function 
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(1.4) 


1  1+^12 


■'^(^1  +  I  £2  -  -S  '^ 


Here   r-,   and  r„   are  position  vectors  of  the  electrons,   b 

Is  the  distance  between  protons,  and  'P-^^     exchanges  r-j^  and 

r^.   We  want  a  function  which  approximates  the  ground  state 
~2 


elgenfunctlon  of  H  almost  as  well  as 
corresponding  elgenfunctlon  of   H^ 
We  Introduce  notation 


approximates  the 


o 


(1.5) 
(1.6) 

(1.7) 


R  -  ^(r.  +  r„) 


r  =  r-,  -  r^ 


^el'-'  =  F  -  ^eo'-'  =  -"'"<F  "  ^> 


We  define  a  basic  set  of  functions  of   r 


(1.8) 


-v2+Vg-L(r) 


f(k^;r)Y.  (0f)  =  k^f(k^;r)Y.  (0|;) 


im' 


Here   (r0^)   are  spherical  coordinates  of   r  ,  Y^  ( 0^)   Is  a 

2 
spherical  harmonic,   V   is  the  Laplaclan  operator  with  respect 

to  coordinates  of   r  ,  and  f(ki;r)  Is  a  regular  solution  of 
(1.8).   To  define  the  nornnallzatlon  of   f  ,  we  use  the  conven- 
tional condition  that  for  values  of  (kr)  much  greater  than 
unity 
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sin[kr  -  i-  ^tt  +  il(ki)] 
(1.9)        f(ki;r)  -  — — ■  •  gy 

Here  •n(k^)   Is  Independent  of  r  and  is  the  usual  phase  shift. 
The  set  f(ki;r)  1.  n    (0^)   is  complete  if  the  operator  on  the 
left  of  (1.8)  has  no  bound  state.   This  holds  for  all  I      if 
it  holds  for  i  -   Q    ,    and  a  sufficient  condition  for  the  non- 
existence of  a  bound  state  is  equivalent  to  a  necessary  condi- 
tion for  the  existence  of  one. 


00 


(1.10)       /   r 
^^o 


-\l(^) 


dr  <  1 


He 


re  r   is  the  point  where  V  -,  (r)   vanishes,  namely   (I/m,). 


(1.11)       -^  exp(-A/M.)  <  1 


Condition  (l.ll)  is  imposed  upon  all  parameter  sets   (A,m.)  .   A 
unitary  transformation  U  is  defined  explicitly 

(1.12)       Uj^(kr)  Y^^(0^)  -  f(k^;r)  ^ ^^^^Vi 

An  operator  equation  follows  from  (1.12) 
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(1.13)       U  ^ 


-V^  +  V  , (r) 
el 


U  =  -V^ 


The  basic  assumption  of  this  paper  Is  that   UO   Is  the  desired 
approximate  elgenfunction  of   H.   This  Is  written  out  ex- 
plicitly using  Fourier  expansions 


(1.14)       *  =  -2  Y.  (0^)  f    A(kim;R)j-  (kr)  k^  dk 

^  Xm  ^"^  -  ^o       ^      ^ 


(1.15)      U*  =  ^  2  Y.  (ef)  f    A(kim;R)f (k^;r)r  k^  dk 

^  im   ^"^  -  ^o 


Here  equation  (l.l4)  simply  defines  the  expansion  coefficients 
A(k^m;R)  ,  and  (1.15)  Is  the  explicit  form  of   U*.   The  analogy 

to  the  original  definition  of  the  Impulse  approximation  Is 

(11) 

clear. 

(1.16)  E  S  =  (4),H^*)  +  r 

(1.17)  S   =  (U4>,U<t>)  =  ((t),4)) 

(i.i8)        r  =  (U*,  (V^  +  Vg^)  U*)  -  (1^,  (V^  +  Vg^)  *) 

In  principle  the  requirement  that  E   Is  a  minimum  determines 
the  parameters  in   *  and  "^g^^^)   simultaneously. 


_  9  _ 


Unfortiinately  the  Integrals  denoted  by  \~'      are  difficult  to 
evaluate.   The  part  Involving  V.   is  expected  to  be  most 

Important  and  its  physical  meaning  is  clear  from  the  impulse 

(12) 
approximation.      The  commutator  of  U  with  V,  ,   which 

determines  the  relevant  part  of   P   ,  is  determined  by  the 
variation  of  V.   over  domains  in  which  V  -,(r)   is  signi- 
ficantly different  from  zero.   If  this  variation  is  large, 
will  be  large  and  positive.   This  is  true  even  if  H  and  H 
have  exactly  the  same  eigenvalues  and  if   *  is  also  the 
exact  eigenfunction  of   H  .   The  point  is  that  in  this  case 
U*  would  be  a  poor  approximation  to  the  correct  eigenfunction 
of  H.   For  this  reason,  the  method  is  not  expected  to  be  use- 
ful for  the  helium  atom,  where  the  maximum  probability  of 
finding  two  electrons  close  is  near  the  origin,  where  V. 
varies  rapidly.   In  the  hydrogen  molecule,  the  corresponding 
probability  has  its  maximum  at  a  saddle  point  of  V „ .   The 

part  of   I   which  Involves  V    is  not  directly  Involved  in 

eo 

limiting  the  range  of  V   ,  and  is  expected  to  be  less  im- 
portant. 

An  approximation  formally  independent  of  the  model,  but 
physically  related,  is  that  the  difference  between  U4'  and 
<J>  involves  only  terms  with  i  =  0  in  (l.l4)  and  (1.15). 
This  implies  the  centrifugal  barrier  masks  the  short-range 
potential  "^qt  (^)  ^°^   higher  values  of  Z.      But  the  evalua- 
tion of   I    remains  difficult  even  with  this  approximation. 
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II.   NON -VARIATIONAL  CONDITION  ON  THE  PARAMETERS  OF   V^-j^lr) 

If   *  and  UO  are  good  approximations  to  the  ground 
state  elgenfunctlons  of  H   and   H  respectively,  the  usual 
manipulation  Involving  hermltlclty  leads  to 


(2.1)  (<f,Vg^U0)  =  (E  -  E^)(4),U<1))  =  0 

Here   E   and  E  ,  the  ground  state  eigenvalues  of  H   and 
H  respectively,  are  approximately  equal.   If  <i>     and  U*  are 
expanded  In  spherical  harmonics  as  in  (l.l4)  and  (l.l5)j  only 
terms  with  £   =  0     are  assumed  to  be  important.   This  is  a 
consequence  of  the  short-range  character  of  "Vg-i(]^)  >    and  is 
equivalent  to  the  approximation  mentioned  in  the  last  para- 
graph of  Section  1.   Direct  substitution  leads  to 

o   ^ 

(2.2)  (-)   /  A(k';R)A(k;R) (k' |V  J  |k)k'^dk'k^dk  d^R  -  0 

IT      J  ~     ~     '  ei 


(2.3)        (k'lVg^llk)  =  /   j^(k'r)V^^(r)f(k;r)r^  dr 

'-'  o 

We  use  the  notation  of  suppressing  the  indices  { £ ,m)      when 
only   ^  =  0   Is  Involved.   The  integral  on  the  right  of  (2.2) 
is  simpler  than  the  integrals   |    defined  by  (l.l8),  but 
further  simplifications  are  needed  to  evaluate  it.   One  of  the 
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simplest  approximations  Is  to  expand   (k'|V  -,||k)   In  powers  of 

(k  -  k  )   and   (k'   -  k  )  ,  keeping  only  powers  zero  and  one, 

2  2 

where   k   Is  some  average  value  of   k  .   This  expansion  Is 

used  In  the  theory  of  scattering  by  short-range  potentials, 

(13) 


such  as  the  "effective  range"  theory, 
becomes 


Condition  (2.2) 


(2.4) 


sin  n(i<^^) 


(k  |V  J  |k  )  = P- 

o  '  el '  '  o         k 

o 


0 


(2.5) 


k2  = 
o 


(-)   r  A(k';R)A(k;R)  k^k'^dk'k^dk  d^R 


(-)   /  A(k';R)A(k;R)  k'^dk'k^dk  d^R 


(2.6) 


0  (R) 

O  'v 


-V^ct)(r,R) 


d^R 


k  = 

o 


4)  (R) 
o  ~ 


d^R 


The  relationship  of  the  phase  shift  to   (k  |V  ^.N^^)   In  (2.4) 
Is  standard/"""^  and  (2.6)  follows  from  (2.5)  If  (l.l4)  and 
the  values  of   j„(kr)   at  r  =  0  are  used.   The  subscript  zero 
on   <I'^(R)   and  the  bracket  in  the  niunerator  of  (2.6)  Indicates 

o  '^ 

that  r  =  0.   We  regard  (2.4)  as  the  definition  of  k   and 

(2.6)  as  Its  evaluation.   We  use  Wang's  values  of  '<  and  b  ; 

(7) 
this  value  of  b   Is  close  to  the  experimental  one. 
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2 
Substituting  (1.4)  Into  (2.6),  we  get  the  value  of   k^ 


„     „  1  +  2Kb  -  ^  (2l<b)^ 

(2.7)        k^  =:  2(<2 ^ 

°       1  +  2iCb  +  4  (2Kb)^ 


We  assume  values  of  A  and  compute  values  of  p.   satisfying 
(2.4)  and  (2.7)  by  the  method  described  in  the  appendix. 
Since   (l/x)   is  essentially  the  range  of  "^g]_(^)  '    the  basic 
theory  requires   (Ab)   somewhat  greater  than  unity  and  this 
guides  the  choice  of  values  of   A.   We  denote  the  Rayleigh-Ritz 
estimate  of  the  eigenvalue  of  H   by  E   ,  and  have 

(2.8)  E^  =   E,  -  (4),V^^4))/S 

(7) 

where   E   Is  Wang^s  estimated  energy^  namely  3«78  ev.     The 
w 

evaluation  of  S  from  (1.4)  and  (1.17)  is  straightforward. 
The  integral   (*,  ^pi  *^   "^^y  ^®  written  as 

(2.9)  (*,  Vg-^^t))  =  (CD,  ^^—   *)  +  n  ^-^(4.,  ^^^_  *) 


The  second  integral  on  the  right  of  (2.9)  is  thus  obtained  from 
the  first  by  differentiation.   The  first  integral  may  be  written 
as  a  sum  of  two  terms,  which  are  counterparts  of  the  usual 
coulomb  and  exchange  integrals  of  molecular  physics,  differing 
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in  that   (e   /r)   replaces   (l/r).   The  direct  part  is  easy 
because   (e   /r)   is  the  Yukawa  kernel.   We  evaluated  the 
exchange  part  by  means  of  the  Gegenbauer  addition  theorem, 
which  is  known  to  quantum  chemists  as  the  Barnett-Coulson 
expansion.   -^   The  usual  method  of  evaluating  the  exchange 
part  by  using  elliptic  coordinates      was  not  used  because 
the  expansion  of  the  Yukawa  kernel  involves  special  functions 
much  more  complicated  then  the  Legendre  functions  which 
occur  in  the  expansion  of  the  Laplace  kernel. 

Computed  values  of  the  second  term  on  the  right  of  (2.8) 
are  shown  in  the  table.   The  variation  with  A   of  tl^is 
integral,  and  with  it  the  expectation  value  of  H   ,  is  some- 
what greater  than  desired.   Since  the  basic  model  seems  good 
intuitively,  we  believe  the  above  approximate  evaluation  of 
k   by  (2.6)  is  responsible  for  the  poor  results  rather  than 
the  use  of   U4)  as  an  Improved  trial  function.   Another  reason 

for  believing  (2.6)  is  a  poor  value  for  k   is  that  (2.7) 

2 
leads  to  negative  values  of   k   for  large   b.   There  is  no 

reason  why  the  basic  model  is  worse  at  large  nuclear  separa- 
tions than  at  equilibrium,  providing  «■      is  replaced  by  unity. 
It  would  be  interesting  to  study   (k'|V  .,  |  |k)   in  detail  and 

see  how  well  it  can  be  approximated  by  a  linear  function  of 

2         2 
k   and  k'    for  all  values  of  these  arj^uments  which  make 

significant  contributions  to  (2.2),  and  for  the  relevant 

domain  of  parameters   (A, p.). 
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APPENDIX 


A  variational  expression  for  m-   Is  derived  from  the 

radial  equation  by  a  straightforward  generalization  of  methods 

(17) 
used  previously.      The  radial  equation  is 


(Al)  ^  ,  (^2  .  V^,)u  =  0 

dr 


where  u  as  usual  is  the  radial  function  multiplied  by  (k  r) . 
Since  the  phase  shift  must  vanish,  u  must  equal  sin  k  r  for 
values  of  r  much  larger  than  the  range  of  V  ,(r).  We  define 
a  function  W  which  vanishes  at  zero  and  Inrinlty 


(A2)         W(r)  =  u(r)  -  sin  k  r 


Direct  substitution  leads  to 


(A3)         dW  ^  (j^2  _  e )^  _  e. gin  k  r  =  -n  e'^^lW+sin  k  r) 

,2or        r        o  o 

dr 


If  W   is  a  function  of  r  which  vanishes  at  zero  and  Infinity, 
is  everywhere  continuous,  and  its  derivative  has  only  a  finite 
number  of  finite  discontinuities,  a  stationary  functional  for  |i. 
is 
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(a4)        [1, 


(^)  _  k^  wf  +  ^ (W,  +  sin  k  r)^ 

dr      o  t    r     t        o 


dr 


00  O     "X 

(W,  +  sin  ^r)      e~        dr 


o 


This  functional  Is  useful  In  a  variation-Iteration  estimate  of 
\i.      We  start  with  an  estimate  of  [x  ,  which  In  practice  Is  the 
Born  approximation,  and  solve  (A^)  numerically.   We  start  one 
solution  at   r  =  0  and  another  at  very  large   r  ,  and  join 
them  continuously  at  the  point   r^  =  (i/a).   We  set  W^   equal 
to  this  function  and  use  (A4)  to  obtain  a  new  estimate  of  \i. 
Then  the  new  value  of  \i      Is  used  In  (A3)  and  the  method  Is 
repeated  until  successive  variational  estimates  of  \i     agree 
to  within  rounding  errors.   In  practice  convergence  Is  very 
fast  and  only  three  or  four  Iterations  are  required.   The 
variational  estimates  of  \i     are  slightly  lower  than  the  Born 
approximation,  and  the  ratio  of  these  two  approximations  of  m- 
Is  usually  0.97  ±  O.OO5.   The  effect  of  using  Born  approxi- 
mation Instead  of  the  variational  estimates  upon  the  expecta- 
tion value  of  V  -,  (r)   Is  larger.   For  example,  for  A  =  1.5 
atomic  units,  the  Born  approximation  leads  to  a  value  of 
0.25  e.v.  for  the  expectation  value  of  V  -,(r)  ,  while  the 
variational  estimate  of  \x      leads  to  O.31  e.v.  ,  as  shown 
In  the  table. 
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Table 

Reduction  In  Zero-Order  Energies  As  a  Function  of 
Range -Determining  Parameter 

A  Expectation  Value  of  V  -.(r) 

atomic  units  electron-volts 

1.0  0.5^ 

1.5  0.31 

2.0  0.20 
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